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i  the  ellipsoid,  thrusting  is  not  necessary.  Variational  performances  in  lateral  thrusting  are  examined  for  a 
hypervelocity,  exoatmospheric,  orbital  vehicle  in  the  final  30  seconds  of  flight  while  it  is  attempting  to 
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indicating  that  the  original  formulation  effectively  uses  the  state  error  estimates. 
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introduction 

CERTAINTY  control'  enhances  interceptor  performance 
by  using  a  terminal  guidance  law  that  incorporates  the 
dynamics  of  the  interceptor  and  target  plus  the  error  knowl¬ 
edge  of  their  estimates.  This  is  done  by  constraining  the  final 
estimated  state  to  a  spherical  inequality  based  on  the  projected 
estimate  error.  The  control  law  reduces  intercept  maneuvering 
when  the  controls  associated  with  cost  do  not  affect  state 
estimate  certainty.2 

Two  variations  of  certainty  control  are  presented  in  an  at¬ 
tempt  to  further  improve  interceptor  performance:  the  first 
uses  a  control  effectiveness  ratio  to  regulate  thrusting  limes; 
the  second  changes  the  certainty  control  constraint  to  an  ellip¬ 
soidal  function  based  on  projected  estimate  error.  Conceptu¬ 
ally,  the  second  variation  produces  a  shrinking  ellipsoid  about 
the  predicted  impact  point  with  the  surface  being  a  function  of 
estimated  error;  if  the  predicted  miss  is  inside  or  touching  the 
ellipsoid,  thrusting  is  not  necessary. 

Variational  performances  in  lateral  thrusting  are  examined 
for  a  hypervelocity,  exoatmospheric,  orbital  vehicle  in  the 
final  30  s  of  flight  while  it  is  attempting  to  intercept  a  boosting 
missile.  System  modeling  and  measurement  processing  are 
identical  to  those  used  in  Ref.  1.  Target  tracking  is  accom¬ 
plished  with  a  ranging  device  and  line-of-sight  sensors  for 
in-plane  and  out-of-plane  measurements.  Noise-corrupted 
data  are  processed  through  an  eight-state  extended  Kalman 
filter  (EKF)  with  serial  updates  occurring  every  0.1  s.  Velocity 
changes  are  determined  by  varying  impact  conditions  using 
splines  to  reduce  computational  burdens  and  allow  a  solution 
that  lends  itself  to  deterministic  techniques. 

Optimum  Spacing  of  Corrective  Thrusts 

Corrective  thrusting  in  me  presence  of  state  estimate  errors 
can  be  optimally  spaced  to  reduce  fuel.  A  control  effectiveness 

Received  Aug.  5,  1991;  revision  received  Oct.  29.  1991;  accepted 
for  publication  Nov.  16.  1991.  This  paper  is  declared  a  work  of  the 
U.S.  Government  and  is  not  subject  to  copyright  protection  in  the 
United  States. 

•Associate  Professor,  Deputy  for  I.abs  and  Research,  Department 
of  Astronautics.  Senior  Member  A1AA 


i  111  IlMNl'J  VOI  !',  NO  4  IM.IM  I  KINt,  NUHS 


ratio,'  p,  is  established  to  determine  the  spacing  between 
thrusts;  the  ratio  directly  yields  thrust  times  when  control  ef¬ 
fectiveness  is  a  linear  function  of  time.  Control  effectiveness  is 
measured  by  the  amount  the  end  conditions  vary  lor  a  speci¬ 
fied  control;  for  this  problem,  as  intercept  time  decreases,  so 
does  the  ability  of  lateral  thrusting  to  vary  miss  distance.  W  ith 
a  ratio  of  two  (p"2.0),  (he  corrective  thrusting  should  only 
occur  when  the  control  has  half  the  effect  ( I /  p)  of  the  pre\  ious 
corrective  thrust.  If  control  effectiveness  is  a  near-linear  func¬ 
tion  of  time,  as  is  the  case  for  a  hypervelocity  orbital  vehicle, 
then  it  will  be  halved  at  about  half  the  time  to  impact  since  the 
last  thrust.  Thrusting  will  occur  at  the  start  of  the  intercept,  at 
one-half  time-to-go,  one-fourth  time-to-go,  one-eighth  time- 
to-go,  etc  When  spacing  is  less  than  the  estimator’s  cycle  time 
(0.1  s  for  this  study),  impact  is  imminent  and  thrust  is  termi¬ 
nated. 


Certainly  Control  Formulation 
with  Ellipsoidal  Constraint 

As  stated  earlier,  if  the  controls  associated  with  cost  do  not 
affect  state  estimate  certainty,  then  fuel  may  be  conserved 
by  using  that  certainty  to  reduce  control  efforts.  When  the 
controls  are  linked  to  the  estimate  certainty,  a  near-perfect 
estimate  yields  the  optimal  control  (certainty  equivalence  solu¬ 
tion'1),  and  a  poor  estimate  causes  a  reduction  in  control.  Cer¬ 
tainty  control  does  this  by  constraining  the  final  estimated 
states  to  a  spherical  inequality  based  on  the  projected  estimate 
error.  An  ellipsoidal  inequality  is  introduced  by  establishing 
the  cost  function. 


AT;  +  A  Vz 
L  =  _ 

subject  to  the  following  constraint: 


(I) 


(i  e.,  attempting  a  /ero-nu\s  solution)  results  111  counterpro¬ 
ductive  maneuvering.  As  the  estimates  improve,  the  sonsuamt 
lightens  and  the  ellipsoid  shrinks,  along  with  predicted  mi’" 
distance.  The  spline  representations  allow  this  stochastiv  pi  oh 
lent  to  be  solved  in  a  deterministic  fashion'  by  adjoining  the 
constraint  to  the  cost  function  to  fotm  the  Hamiltonian 
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where  AT„  and  AT.  are  the  interceptor's  velocity  changes,  and 
K  is  the  constraint  weighting  factor.  Time-to-go  rgo  is  used 
as  a  third  control  parameter  to  minimize  miss  distance  but  does 
not  explicitly  appear  in  the  cost  function.  The  final  state  esti¬ 
mates  and  their  deviations  (oxf,oV),o-j)  are  deter¬ 

mined  by  running  the  filter  forward  to  predicted  impact  time 
without  measurement  or  control  updates  and  then  representing 
their  time  history  with  splines; 


Equations  (2),  (13),  (14),  and  (15)  constitute  lour  equations 
with  four  unknowns,  which  can  be  reduced  to  two  equations 
and  two  unknowns  using  Eqs.  (7)  and  (8).  Substituting  Eq.  (7) 
into  Eq.  (13)  yields 
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■>  . 

In  a  similar  manner,  substituting  Eq.  (8)  into  Eq.  (14)  yields 


x>  -  Axtga  +  Bxt‘0  +  Cxt g0  +  D,  (3) 

Ts  =  A  ,./g0  +  Bv/g0  +  Cytto  +  Dv  (4) 

*5  =  A  4  'go  +  B.  /g20  +  C,  /go  +  D.  (5) 
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Equations  (2)  and  (15)  can  now  be  solved  in  terms  of  X  and 
/go,  with  AT,  and  AT.  determined  afterward  from  Eqs.  (22)  and 
(24).  The  parameters  X  and  tf„  can  be  found  by  numerical 
techniques  using  the  Jacobian; 
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Conceptually,  the  constraint  produces  a  deviation  ellipsoid 
about  the  predicted  impact  point.  If  the  predicted  miss  is  inside 
or  touching  the  ellipsoid,  then  thrusting  is  not  necessary.  If  the 
predicted  miss  is  outside  the  ellipsoid,  then  minimum  thrusting 
is  determined  to  bring  the  miss  to  the  ellipsoid’s  surface.  Sen¬ 
sor  inaccuracies  will  cause  the  predicted  impact  point  (ellip¬ 
soid  center)  to  jitter  with  each  measurement;  chasing  this  point 
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If  the  stales  are  perfectly  known,  the  a  terms  will  he  zero 
and  the  constraint  of  E-.q.  (2)  will  only  be  satisfied  with  a  ' 
predicted  miss  of  zero;  the  ellipsoidal  certainty  control  equa¬ 
tions  reduce  to  the  certainty  equivalence  optimal-control  for¬ 
mulation.1  If  the  estimate  is  poor,  the  o  terms  will  be  la.ge  and 
the  inequality  constraint  of  Eq.  (2)  will  result  in  very  little,  if 
any,  change  in  velocity. 

Computer  Simulation 

A  head-on,  10-deg  out-of-plane  intercept  is  examined  w  ith  a 
time-to-go  of  30  s.  The  interceptor  is  initially  traveling  at  12 
km/s  at  an  altitude  of  750  km  with  a  lateral  acceleration  range 
of  3-60  m/s:  in  each  axis.  The  booster’s  initial  acceleration  is 
3.15788  m/s:,  with  a  unitized  mass  flow  rate  of  0.01579s  '. 

A  time  lag  of  0. 1  s  is  used  when  computing  velocity  changes 
to  account  for  measurement  processing,  controller  process¬ 
ing,  and  thruster  response.  Target  acquisition  is  assumed  to 
take  3  s;  thrusting  is  not  permitted  during  this  time.  This  sim¬ 
ulation,  written  in  Fortran  77  to  run  on  a  VAX  3600,  generates 
200  Monte  Carlo  runs  per  case. 
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Fig.  1  Performance  of  certainly  control  for  head-on,  10-deg  out-of¬ 
plane  intercept. 
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Fig.  2  Performance  of  certainty  control  for  head-on,  10-deg  out-of- 
plane  intercept  with  p  =  1.1. 
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10-deg  out-of-plane  intercept. 


Results 

The  original  certainty  control  performance  serves  as  a  basis 
of  comparison  for  the  variations,  with  Fig.  1  showing  the 
relationships  of  miss  distance  and  total  AF  to  the  constraint 
parameter  K .  Figure  2  illustrates  the  relationships  when  the 
thrusting  times  are  limited  by  a  control  effectiveness  ratio  of 
1.1,  and  Fig.  3  shows  the  effect  of  altering  the  control  con¬ 
straint  from  spherical  to  ellipsoidal. 

The  minimum  AF  performances  of  all  three  control  strate¬ 
gies  are  found  in  Table  1 .  A  control  effectiveness  ratio  of  1,1 
slightly  reduces  the  miss  distance  and  increases  the  total  AT, 
whereas  larger  ratios  degrade  performance  for  this  particular 
intercept.  Reformulating  the  control  law  using  an  ellipsoidal 
constraint  reduces  the  total  AF,  but  sacrifices  some  accuracy. 
Although  they  are  inconclusive,  these  results  indicate  that  the 
original  formulation  of  certainty  control  effectively  uses  the 
stale  deviations  to  minimize  maneuvering  costs  while  main¬ 
taining  a  high  level  of  accuracy. 

Conclusions 

In  this  Note,  two  variations  of  certainty  control  were  exam¬ 
ined  to  determine  their  c-^nbilKv  to  mir'mize  iatcr„l  .eloci'y 
changes  ol  a  iiypervelocuy  orbital  veiiicle  in  a  head-on.  10-deg 
out-of-piane  intercept.  The  first  variation  used  a  control  effec¬ 
tiveness  ratio  to  regulate  thrusting  times;  the  second  changed 
the  spherical  constraint  function  to  an  ellipsoidal  one.  Neither 
variation  simultaneously  reduced  maneuvering  cost  and  miss 
distance  when  compared  to  the  original  formulation  of  cer¬ 
tainty  control,  indicating  that  the  original  formulation  of 
certainty  control  effectively  uses  the  state  deviations  to  mini 
mize  maneuvering  costs  while  maintaining  a  high  level  of  accu¬ 
racy.  The  ellipsoidal  variation  best  demonstrated  a  tradeoff  in 
accuracy  to  reduce  maneuvering  costs,  a  choice  to  be  made 
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based  on  mission  constraints.  An  area  ot  further  research  is  to 
find  a  way  to  make  better  use  of  the  stale  deviations  in  Vor mu¬ 
tating  the  control  law,  i.e.,  reducing  cost’ while  further  improv¬ 
ing  accuracy. 
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